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Abstract. We report the discovery of a massive planetary companion orbiting the young disc star HD 81040. Based on five 
years of precise radial-velocity measurements with the HIRES and ELODIE spectrographs, we derive a spectroscopic orbit 
with a period P = 1001.0 days and eccentricity e = 0.53. The inferred minimum mass for the companion of m2 sini = 6.86 
Mjup places it in the high-mass tail of the extrasolar planet mass distribution. The radial-velocity residuals exhibit a scatter 
significantly larger than the typical internal measurement precision of the instruments. Based on an analysis of the Ca II H and 
K line cores, this is interpreted as an activity-induced phenomenon. However, we find no evidence for the period and magnitude 
of the radial- velocity variations to be caused by stellar surface activity. The observed orbital motion of HD 81040 is thus best 
explained with the presence of a massive giant planet companion. 
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1. Introduction 

Except for one experiment (Cochran et al. 2002), radial- 
velocity surveys for planets orbiting nearby, solar-type stars 
with a targeted precision of 1 - 5 m s"' (e.g.. Mayor et al. 
2004; Santos et al. 2004; Butler et al. 2004; Marcy et al. 
2005a, 2005b) have systematically attempted to exclude from 
their samples chromospherically active stars. This is because 
changes in the visibility of active regions and variations of 
the steUar absorption line profiles can possibly cause signif- 
icant deterioration in the achievable radial-velocity precision. 
However, for about a dozen planet hosts significant activity 
levels have not prevented radial-velocity surveys from detect- 
ing Keplerian signals due to orbiting giant planets. We report 
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in this paper on radial-velocity measurements of the young disc 
star HD 81040 (HIP 46076, BD-l-20°2314). These observations 
reveal the presence of a massive planet candidate orbiting the 
star, with a minimum mass of 7.61 Mjup (where Mjup is the 
mass of Jupiter). 

The variable velocity of HD 81040 was first detected by the 
G-Dwarf Planet Search (Latham 2000), a program designed to 
conduct a first reconnaissance for giant planets orbiting a sam- 
ple of nearly 1000 nearby G dwarfs. This survey employed the 
HIRES spectrograph (Vogt et al. 1994) on the 10-m Keck 1 
telescope at the W. M. Keck Observatory (Hawaii). Follow- 
up observations were carried out within the context of the 
ELODIE Planet Search Survey (Mayor & Queloz 1995; Perrier 
et al. 2003), which uses the ELODIE fiber-fed echelle spec- 
trograph (Baranne et al. 1996) on the 1.93-m telescope at the 
Observatoire de Haute-Provence (CNRS, France). 

Radial velocities with ELODIE are obtained by cross- 
correlating the observed spectra with a numerical template. 
The "simultaneous thorium-argon technique" with dual fibers 
(Baranne et al. 1996) allows for the monitoring and correc- 
tion of instrumental drifts. The precision achieved with this 
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instrument on bright, inactive stars is ~ 6.5 m s ' (Perrier et 
al. 2003). The HIRES instrumental profile and drifts are moni- 
tored using an iodine absorption cell (Marcy & Butler 1992). In 
order to derive radial velocities, the data reduction procedure 
involves the modelling of the temporal and spatial variations 
of the instrumental profile of the spectrograph (Valenti et al. 
1995), and is conceptually similar to that described by Butler 
et al. (1996). Internal errors, computed from the scatter of the 
velocities from the echelle orders containing I2 lines, are typi- 
cally 8-10 m s"' for inactive, solar-type stars, for the relatively 
short exposures adopted for the G-dwarf Planet Search. 

The collaboration between our two teams has so far resulted 
in three important discoveries: the determination of the spectro- 
scopic orbit of the first extrasolar planet transiting the disc of 
its parent star (HD 209458b; Mazeh et al. 2000), the detection 
of HD 80606b, the planet with the largest eccentricity known 
to-date (e = 0.927; Naef et al. 2001a), and the identification of 
the first giant planet in a stellar triple system (HD 17891 IBb; 
Zucker et al. 2002). 

The Keck/HIRES observations of HD 8 1040 began in April 
1999 and have a time baseline of ~ 8.5 months. The ELODIE 
follow-up observations span ~ 3.3 years, starting in February 
2002. We present in Sect. 2 the radial-velocity data and the 
orbital solution derived utilizing the combined datasets. The 
stellar characteristics of HD 81040 are considered in Sect. 3. 
Finally, Sect. 4 is devoted to a summary and discussion of our 
findings. 

2. Radial-Velocity Data and Orbital Solution 

The 26 radial velocity measurements for HD 81040 are listed 
in Tabled In order to bring the HIRES data to the ELODIE 
system, they have been initially shifted by an amount corre- 
sponding to the J velocity measured by ELODIE, ISRY - 49.24 
km s"^ A residual velocity offset ARVh-e between the two 
is included as an additional free parameter in the best-fit or- 
bital solution presented in Table |2] Formal uncertainties in the 
parameters are obtained by re-scaling by the reduced x'' of 
the solution. We show in the upper and lower panel of Fig. 
the radial-velocity measurements as a function of time and or- 
bital phase, respectively. Post-fit velocity residuals are shown 
in the two sub-panels. Using the value for the primary mass of 
M^^ = 0.96 + 0.04 Mq obtained as described in Sect. 3, we de- 
rive a minimum mass for the companion m2 sin / = 6.86 + 0.71 
Mjup. The radial-velocity variations measured in HD 81040 
thus reveal the presence of a long -period (P ^ 2.7 years), ec- 
centric (e = 0.53), massive planetary-mass object. Due to gaps 
in phase coverage and limited number and time baseline of the 
observations, further measurements will be required to improve 
on the determination of some of the orbital elements (such as 
the orbital period and the eccentricity). 

Given the values of sin ; and P, translating into an as- 
trometric signature a ~ //as at the distance of HD 81040, 
we have investigated the Hipparcos Intermediate Astrometric 
Data (IAD) in an attempt to further constrain the companion 
mass (for a review of the method see for example Sozzetti 
2005, and references therein). The star is flagged as single in 
the Hipparcos database. The best-fit astrometric solution with 



Table 1. Radial velocity measurements of HD 81040. The 
HIRES (H) velocities are shifted to the ELODIE (E) zero point. 



BJD - 2,400,000 


Radial Velocity 


0"RV 


Instrument 




(km s-') 


(kms-') 




51291.36542 


49.319 


0.012 


H 


51293.37376 


49.312 


0.014 


H 


51545.63921 


49.088 


0.009 


H 


52308.63310 


49.399 


0.012 


E 


52356.41490 


49.412 


0.021 


E 


52359.45970 


49.363 


0.012 


E 


52360.42310 


49.391 


0.011 


E 


52615.71400 


49.119 


0.011 


E 


52649.66360 


49.126 


0.011 


E 


52649.67630 


49.129 


0.013 


E 


52719.40970 


49.100 


0.008 


E 


52722.49520 


49.118 


0.012 


E 


52993.63610 


49.204 


0.012 


E 


53034.56760 


49.293 


0.015 


E 


53094.36570 


49.308 


0.010 


E 


53101.39950 


49.285 


0.010 


E 


53359.69070 


49.389 


0.016 


E 


53361.66000 


49.412 


0.021 


E 


53421.52680 


49.459 


0.019 


E 


53428.48940 


49.462 


0.026 


E 


53461.36050 


49.414 


0.016 


E 


53463.41260 


49.458 


0.020 


E 


53464.39850 


49.408 


0.011 


E 


53486.35380 


49.359 


0.011 


E 


53491.32410 


49.348 


0.010 


E 


53518.35250 


49.234 


0.013 


E 



Table 2. Best-fit orbital solution for HD 81040, derived mini- 
mum companion mass, and properties of the post-fit residuals. 



Parameter 


Value 


P (days) 


1001.7 ±7.0 


T (BJD) 


2452504.0 ± 12.0 


e 


0.526 ± 0.042 


y (km s"') 


+49.2535 ± 0.0063 


ojn 


81.3 ±7.2 


Ki (m s-') 


168 ±9 


ARVn-E (kms-') 


+0.057 ± 0.019 


ai sini (xlO* km) 


1.96 ±0.19 


Mm) (xlO-'Me) 


3.0 ±0.9 


m2 sin / (Mjup) 


6.86 ±0.71 


a(AU) 


1.94 


A' 


23(E) + 3(H) 


< ctrv > (tn s"') 


13.7 


a-Q-c (ms"') 


26.0 (E: 25.0, H: 6.0) 




3.26 


Prix^) 


2 X 10-' 



the period of the spectroscopic orbit gives a statistically in- 
significant semi-major axis of 2.0 +1.0 mas. We have then 
used these values to derive an upper limit to the actual size of 
the astrometric orbit a with a 2.3cr confidence, following the 
procedure described in Zucker & Mazeh (2001). The resulting 
upper limit on the companion mass of ~ 75 Mjup, at the 99% 
confidence level, clearly indicates its sub-stellar nature. 
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Fig. 1. Top: Radial-velocity measurements as a function of time 
for HD 81040. Open circles identify Keck/HIRES velocities, 
while filled black circles coiTespond to the ELODIE dataset. 
Bottom: radial velocities as a function of orbital phase. The 
two sub-panels show the post-fit residuals in both cases. 

The scatter o-q-c of the post-fit residuals reported in Tabled 
is abnormally large when compared to the average < ctrv > of 
the internal errors for both instruments. Inspection of the veloc- 
ity residuals (top and bottom sub-panels of Fig. ^ by means of 
a periodogram search revealed no additional periodicity within 
the time-span of the observations. No significant velocity trend 
can also be found in the post-fit residuals (a linear fit to the 
residuals time series had an insignificant slope of 1.8 + 2.2 
m s ' yr '). Thus, another companion on a longer-period or- 
bit superposing a second radial-velocity signal is an unlikely 
explanation (although further measurements will help to better 
address this question). Another possibility to explain the excess 
jitter in the residuals is to invoke effects due to chromospheric 
activity of the parent star (e.g., Baliunas et al. 1995; Saar et al. 
1998; Santos et al. 2000a). We discuss this possibility in the 
next Section. 



3. Properties of the Host Star 

HD 81040 (HIP 46076, BD-h20°2314) is a bright, nearby 
dwarf. Hipparcos astrometry places it at 32.56 pc from the 
Sun. The quoted visual magnitude and color index from 
Hipparcos are my - 7.72 and B - V - 0.68, respectively. 



Table 3. Observed and infeiTed stellar parameters for HD 
81040. Spectral type, apparent magnitude, colour index, paral- 
lax, and proper motions are from Hipparcos (ESA 1997). The 
atmospheric parameters Tetj, \ogg, ^t, and [Fe/H], and the stel- 
lar mass and radius have been derived as described in Sect. 3.1. 
The values of v sin /, < log R'^^^ >, log e(Li), and age have been 
obtained as described in Sect. 3.2. 



Parameter 


Value 


5 p. Type 


G2/G3 


my 


7.72 


B-V 


0.68 


n (mas) 


30.71 ± 1.24 


d (pc) 


32.56 ± 1.31 


Pa (masyr"') 


-151.35 ± 1.08 


Ps (masyr"') 


35.91 ±0.52 


Mv 


5.17 


B.C. 


0.10 


{U,V,W) (kms-') 


(40.9,-1.5,25.1) 


Teff (K) 


5700 ± 50 


logg (cgs) 


4.5+0.1 


(km s-') 


0.95 + 0.05 


[Fe/H] 


-0.16 ±0.06 


(Mo) 


0.96 ± 0.04 


R* (Ro) 


0.86 ± 0.04 


< logiJ^j, > 


-4.48 


V sin i 


2± 1 


loge(Li) 


1.90/1.91 ±0.07 


< t > (Gyr) 


0.73 ±0.1 



The corresponding absolute magnitude is My - 5.17, con- 
sistent with a G2/G3 spectral type (Cox 2000), and this clas- 
sification is also favored by our effective temperature deter- 
mination (see Sect. 3.1). A GO spectral type is reported in 
SIMBAD and in several works in the literature, however the 
ultimate source is the Henry Draper catalog, in which the spec- 
tral classification is rather coarse. Using proper motions and 
parallax from Hipparcos and the y velocity obtained from the 
orbital solution described in Sect. 2 (in good agreement with 
the Nordstrom et al. (2004) reported value of 48.9 ± 0.2 km 
s '), the Galactic velocity vector for HD 81040 (coiTected for 
the Local Standard of Rest following Mihalas & Binney 1981) 
is then (U,V,W) = (40.9, -1.5, 25.1) km s"'. These values are 
summarized in Tabled 



3.1. Abundance Analysis 

We have utilized the Keck/HIRES high-signal-to-noise tem- 
plate spectrum of HD 8 1040 to carry out an LTE spectroscopic 
iron abundance analysis. In this study, we have followed the 
approach described in previous works (Sozzetti et al. 2004, 
2005b, and references therein). Our final iron list consisted of 
28 Fe I and 6 Fe II lines (taken from Sozzetti et al. 2005b). 
The final atmospheric parameters, summarized in Tabled are: 
Teff = 5700 + 50 K, log^ = 4.5 + 0.1 (cgs), = 0.95 + 0.05 
km s and [Fe/H] - -0.16 + 0.06. All these numbers are in 
very good agreement with those reported by Allende Prieto & 
Lambert (1999) and Nordstrom et al. (2004). 
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The spectroscopically determined values of T^.^ and [Fe/H], 
and the absolute luminosity estimate, were then used as input 
to the Yale stellar evolution models (Yi et al. 2003) to derive 
estimates of the stellar mass and radius for HD 81040, as well 
as their uncertainties. The results, also reported in Table|3] are 
as follows: = 0.96 ± 0.04 Mo, R^, = 0.86 + 0.04 Rq. The 
predicted log g value is 4.52, in excellent agreement with the 
spectroscopic estimate. 



3.2. Youth and Activity Indicators 

Montes et al. (2001), on the basis of its galactic kinematics, 
classify HD 81040 as a young Galactic disc star, with no clear 
membership to any stellar kinematic group. The star had been 
selected in that study based on its chromospheric activity levels 
as measured by Strassmeier et al. (2000) in the Vienna-KPNO 
search for Doppler-imaging candidate stars. For HD 81040, 
Strassmeier et al. (2000) measured a value of the chromo- 
spheric emission ratio logR'^^y^ = -4.52. All our Keck/HIRES 
spectra show significant core reversal of the Ca II H and K 
lines. In the top panel of Fig.|2]we show a region of the HIRES 
template spectrum centered on the Ca II H line. The emission 
feature is clearly visible. For comparison, in the bottom panel 
we show the spectrum of an old, inactive star (HIP 105888) 
with the same temperature from the Sozzetti et al. (2005a) 
sample of metal-poor stars, which has a metallicity of [Fe/H] 
= -0.72. 

We have measured the Mount Wilson chromospheric activ- 
ity index S (Duncan et al. 1991) from the Ca II H and K lines 
in the Keck/HIRES spectra, and converted it to ^^j^, following 
the approach described in Sozzetti et al. (2004). Based on four 
data-points, the average activity level measured in HD 81040, 
< logR'^^ >- -4.48, agrees very well with that reported by 
Strassmeier et al. (2000). Based on the Noyes et al. (1984) em- 
pirical calibrations, the mean inferred chromospheric age for 
HD 81040 is then < t >= 0.73 + 0.1 Gyr, and the average 
stellar rotation period < Prot >= 9.8 days. On the other hand, 
the stellar projected rotational velocity appears low. Its mea- 
sured value, using the mean ELODIE cross-correlation func- 
tion (CCF) dip width, is v sin/ - 2+1 km s in excellent 
agreement with the one (2.0 km s"') reported by Nordstrom et 
al. (2004). Our estimate of v sin / is also compatible with the 
two values of 3.7 and 4.9 km s"' (with typical uncertainties of 
2-4 km s"') quoted by Strassmeier et al. (2000). 

Another important piece of circumstantial evidence in fa- 
vor of a young age for HD 81040 is the presence of a signifi- 
cant Lithium (Li) feature. In order to derive an estimate of the 
Li abundance, we summed all the ELODIE spectra in the re- 
gion of the A = 6707.8 A line. We then utilized the resulting 
co-added ELODIE spectrum to carry out a spectral synthesis of 
the Li line, using the atmospheric parameters derived from the 
Fe-line analysis and the same line list of Reddy et al. (2002). In 
Fig.|3]we show the comparison of the spectrum of HD 81040 
with three different models, differing only in the Li abundance. 
The best-fit model gives as a result log e(Li) = 1.90. Following 
Naef et al. (2001b) and Perrier et al. (2003), we also directly 
measured the equivalent width EW^ of the A - 6707.8 A Li I 
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Fig. 2. Comparison between the Ca II H line for HD 8 1040 (up- 
per panel) and an inactive star of the same temperature (lower 
panel). 
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Fig. 3. A portion of the ELODIE co-added spectrum of HD 
81040 containing the Li I line at 6707.8 A (asterisks), com- 
pared to three syntheses (lines of different colours and styles). 
A significant feature is clearly detected. 

line, and used the Soderblom et al. (1993) curves of growth to 
obtain a second, independent assessment of the Li abundance. 
The resulting value is log e(Li) = 1.91+ 0.07, with the reported 
error having been estimated by changing EW^ and T^s by + 1 cr. 
Thus, the two methods we applied for the determination of the 
Li abundance in HD 81040 are in very good agreement with 
each other, and our numbers are also compatible with the value 
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Fig. 4. Left: Bisector inverse slope of the ELODIE CCFs (as defined in Queloz et al. 2001) as a function of the radial velocities 
obtained from the CCFs. A linear fit to the data (solid line) has a statistically insignificant slope of -0.010(±0.026)yf. A rank- 
correlation test gave a probability of no-correlation p - 0.41. Right: Bisector span as a function of orbital phase (p. A linear 
fit to the data (solid line) has a statistically insignificant slope of -O.OOO9(+OO9)0. A rank-correlation test gave a probability of 
no-correlation p - 0.58. 



of log e(Li) = 2.13 + 0.15 quoted by Strassmeier et al. (2000). 
The estimated Li abundance is consistent with the typical val- 
ues for a star slightly older than the Hyades cluster (in good 
agreement with the chromospheric age estimate) and with the 
same temperature of HD 8 1040, both on an empirical as well as 
theoretical basis (see for example Sestito and Randich 2005). 

3.3. Study of Radial-Velocity Jitter 



extra-scatter cr, 



able cause for the observed radial-velocity curve. To put this 
conclusion on firmer grounds, we have carried out an anal- 
ysis of the line bisectors from the ELODIE CCFs, searching 
for possible correlations between the bisector span and radial- 
velocity measurements, orbital phase, and velocity residuals. 
Rank-correlation tests gave probabilities of no correlation rang- 
ing between 0.41 and 0.92. As an example, in the two panels of 
Fig.|4]we show the bisector span values plotted against the ve- 
locities obtained from the ELODIE CCFs and the orbital phase, 
respectively. No trend is visible, further supporting the expla- 
nation of true stellar reflex motion induced by an orbiting com- 
panion as opposed to radial-velocity variations originating in 



With all the information gathered on the parent star pointing 
in the direction of a young age for HD 81040, the observed 

^|^o-c~ ^ ''"Rv >^ ~ 22 m s 'in the veloc- jjjg stellar atmosphere. Unfortunately, the R'.^^ measurements 

are too few to undertake a similar analysis. Finally, Hipparcos 
lists HD 81040 as photometrically stable, with a scatter of 9 
mmag. No significant signal at any frequency was found as a 
results of a periodogram analysis on the photometric data. 

4. Summary and Discussion 

We have presented in this work a combined Keck/HIRES -I- 
ELODIE dataset of radial-velocity measurements for the young 
disc star HD 81040. The measurements reveal the presence of 
a massive (m2 sin / = 6.86 Mjup) planetary companion, on an 
eccentric (e - 0.53), relatively long-period {P - 1001 days) 
orbit. The host star is a G2/G3 dwarf with a metallicity [Fe/H] 
= -0.16, very close to the average ([Fe/H] ^ -0.1) of the solar 
neighbourhood (Nordstrom et al. 2004), but somewhat metal- 
deficient with respect to the average ([Fe/H] ^ 0.14) of the 
metallicity distribution of planet-hosting stars (e.g., Fischer & 
Valenti 2005). Both the chromospheric activity level (measured 
using the Ca II H and K fines in the HIRES spectra) as well as 



ity residuals can then be interpreted in terms of activity-related 
processes. It is well known (e.g., Saar et al. 1998; Santos et al. 
2000a) that spectral line profile variations induced by surface 
activity (e.g., spots) can translate into excess radial- velocity jit- 
ter. For an early-G dwarf with the level of activity exhibited by 
HD 81040, empirical estimates (e.g., Saaret al. 1998; Santos et 
al. 2000a; Paulson et al. 2002; Wright 2005) indicate a typical 
value of the activity-induced radial-velocity jitter of cr[, ^ 20 m 
s"', in good agreement with the observed scatter in the post-fit 
velocity residuals for HD 81040. 

However, the magnitude and period of the detected Doppler 
signature are such that its interpretation in terms of the pres- 
ence of a massive planetary companion is convincing. In fact, 
the estimated rotation period is about two orders of magnitude 
shorter than the inferred orbital period, and the radial-velocity 
semi-amplitude is almost an order of magnitude larger than the 
presumed activity-induced level of jitter. Unlike the case, for 
example, of the short-period variable HD 166435 (Queloz et 
al. 2001), stellar surface activity alone is thus a very improb- 
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Table 4. Main properties of active planet-hosting stars and their orbiting planets. Only objects with a chromospheric age 
estimate f < 1.5 Gyr are included (corresponding to logR'^^ > -4.65). Columns 1 through 12 report: star name, spectral type, 
metallicity, projected rotational velocity, chromospheric emission ratio, lithium abundance or upper limit (where applicable), 
estimated rotation period, chromospheric age, radial velocity jitter estimate, orbital periods of the companions, minimum 
masses, and radial-velocity semi-amplitude. Unless otherwise noted, all parameter values are taken from the literature sources 
listed in column 12. The star HD 166435 does not have a planet, but is listed for comparison. 



Star 


Sp. 


[Fe/H] 


V sin i 


log RUi, 

HK 


log eCLi)" 




Age'' 


<^ 

V 


P 


1112 sin ; 




Refs."" 




Type 




(m s-') 






(days) 


(Gyr) 


(m s-') 


(days) 


(Mj„p) 


(ms-') 




HD 81040 


G2/G3 


-0.16 


2.0 


-4.48 


1.90 


10 


0.7 


7/19 


1001 


6.86 


168 


(1) 


HD 1237 


G6V 


+0.10 


5.5 


-4.44^' 


2.24 


10 


0.6 


21/21 


133 


3.32 


164 


(2) 


HD 17051 


GOV 


+0.03 


6.1 


-4.65 


2.63 


12 


1.5 


13/14 


320 


2.26 


67 


(3) 


HD 22049 


K2V 


-0.10 


2.5« 


-4.47'' 


< 0.3 


10 


0.7 


9/9 


2502 


0.86 


18 


(4) 


HD 40979 


F8 


+0.19 


7.4 


-4.63 


2.79 


12 


1.0 


19/15 


263 


3.28 


100 


(5) 


HD 73256 


G8/K0 


+0.29 


3.2 


-4.49 




12 


0.8 


13/19 


250 


1.87 


265 


(6) 


HD 121504 


G2V 


+0.16 


2.6 


-4.57 


2.66 


9 


1.1 


7/16 


63 


1.22 


45 


(7) 


HD 128311 


KO 


+0.08 


5.7 


-4.44 


< -0.4 


10 


0.5 


11/9 


458 
928 


2.18 
3.21 


67 
76 


(8) 


HD 130322 


KO 


-0.02 


1.9 


-4.39 


< 0.2 


8 


0.3 


7/9 


10 


1.02 


115 


(9) 


HD 141937 


G2/G3V 


+0.01 


2.1 


-4.65 


2.48 


13 


1.5 


6/14 


653 


9.70 


247 


(10) 


HD 142415 


GIV 


+0.21 


3.3 


-4.55 




10 


1.1 


8/16 


386 


1.62 


51 


(7) 


HD 147513 


G3/5V 


+0.06 


1.5 


-4.38 


2.05 


5 


0.3 


4/23 


528 


1.21 


31 


(7) 


HD 192263 


K2V 


-0.14 


1.8 


-4.39 


< -0.3 


9' 


0.3 


6/9 


24 


0.72 


62 


(11) 


HD 196050 


G3V 


+0.22 


3.1 


-4.65' 


2.15 


16 


1.5 


9/14 


1321 


3.02 


49 


(7) 


HD 166435 


GO 


-0.07 


7.6 


-4.26 


< 1.7 


4" 


0.2 


33/28 


4 


0.60 


83 


(12) 



" All abundances and upper limits are taken from Israelian et al. (2004), except for HD 81040 (this work), HD 40979 (Fischer et al. 2002) 
and HD 166435 (Queloz et al. 2001). 

* From the Noyes et al. (1984) calibrations. 
'' From the Noyes et al. (1984) calibrations. 

The two values are derived using the Saar et al. (1998) and the Santos et al. (2000a) calibrations, respectively. 

" (1): this work; (2); Naef et al. 2001b; (3) Kurster et al. (2000); (4); Hatzes et al. 2000; (5): Fischer et al. 2002; (6): Udry et al. 2003; (7): 
Mayor et al. 2004; (8): Vogt et al. 2005; (9): Udry et al. 2000; (10): Udry et al. 2002; (11): Santos et al. 2000b; (12): Queloz et al. 2001. 

^ From Henry et al. (1996). Naef et al. (2001b) report a substantially higher value log^Jj^. = -4.27 (but with large scatter), corresponding 
to Py„t = 4 days and t = 0.02 Gyr. 

* From Tokovinin (1992). Uncertainties on the v sin / value for HD 22049 are rather large. 
'' From Henry et al. (1996). 

' A photometric period of 24 days, matching the orbital period, was obtained by Henry et al. (2002). However, Santos et al. (2003) have 
shown the planet hypothesis still holds. 
' Henry et al. (1996) report a much lower value logR'j^^ = -5.04. 

* The photometric period matches the orbital and chromospheric period estimates. Queloz et al. (2001) have shown the observed signal 
is due to activity-related processes. 



the Lithium abundance (estimated using the co-added ELODIE 
spectrum) speak in favor of a young age for HD 81040 (~ 0.8 
Gyr), and our results agree well with previous studies of this 
object (Strassmeier et al. 2000; Monies et al. 2001). 

The excess scatter in the residuals from the best-fit orbit 
with respect to the nominal internal errors of the HIRES and 
ELODIE measurements, ~ 22 m s"', is likely attributable to 
the high level of activity of HD 81040. Radial-velocity surveys 
of nearby F-G-K dwarfs tend to avoid chromospherically ac- 
tive stars, for which intrinsic velocity jitter caused by surface 
inhomogeneities (e.g., spots, convection) not only degrades the 
achievable measurement precision, but, when correlated with 
the stellar rotation period, can even mimic the signal produced 
by an orbiting companion (e.g., Queloz et al. 2001). However, 
HD 81040 is not the first young star to have been found har- 
boring a planetary-mass companion, despite its significant ac- 
tivity levels. We summarize in Table |3 the properties of active 



(log/?^j^ > -4.65, corresponding to a chromospheric age esti- 
mate f < 1 .5 Gyr) stars with known detected planets, and the 
main characteristics of the latter. From inspection of the Table, 
at least two important considerations can be made. 

First, the two main youth indicators (logR'^^ and, when 
measurable, log e(Li)), and to a lesser extent v sin /, are in fair 
agreement with each other to indicate relatively young ages for 
these stars. We show in the top and middle panels of Fig. |5] 
a comparison between the average values of v sin / and of the 
chromospheric flux of the Ca II K line core F'^ (expressed in 
erg cm"^ s"') for the ensemble of active planet hosts in Tabled] 
and the age-activity-rotation relationships for open clusters and 
average observed values from Pace & Pasquini (2004), and ref- 
erences therein. The agreement is broad, and a similar result 
holds for the comparison between average Li abundances in 
open clusters from Sestito & Randich (2005) and the average 
log e(Li) for the ensemble of active stars with planets (lower 



A. Sozzetti et al.: A Massive Planet to the Young Disc Star HD 81040 



7 



3.0 - 




n.o h I I I 

1 a 3 4 

Age (Gyr) 



14 p 
12 

10 - 




n h 

12 3 4 

Age (Gyr) 



3.5r 




1.5 - 



O.n h , , , , I , , , , I , , , , I , , , , I , , , , I , , , , I , , , , - 
6.5 7.0 7.5 8.0 8.5 9.0 9.5 10.0 

Log(age) (yr) 

Fig. 5. Top: average chromospheric flux F'^ (corrected for the 
photospheric-flux contribution) as a function of age for four 
intermediate-age clusters (diamonds), the sun (asterisk), and 
active stars with planets (filled circle). For clusters, error bars 
represent peak-to-peak spreads. For stars with planets, they cor- 
respond to the spread around the mean value. Lines of differ- 
ent shapes represent theoretical power-law fits to the data with 
different exponents (Pace & Pasquini 2004). Center: the same, 
but in the v sin /-age plane. Bottom: average Lithium abundance 
log e(Li) as a function of age for open clusters (diamonds) and 
active planet hosts (filled circle). For the cluster averages, only 
G-type stars are considered (Sestito and Randich 2005). For 
stars with planets, only actual log e(Li) estimates (no upper lim- 
its) have been used. All error bars correspond to the dispersion 
around the mean value. 



panel of Fig. |5}- The average age inferred from the chromo- 
spheric activity values for the stars in Table |3 is ~ 0.8 Gyr, 
slightly older than the Hyades cluster, and indeed the average 
values of F'^, v sin /, and log e(Li) are all consistent with such 
interpretation. 



Second, to appreciate how severe a danger (moderate) 
youth and activity can pose for radial-velocity planet surveys, 
it is useful to compare the values of radial-velocity semi- 
amplitude and orbital period of the companions with the ex- 
pected levels of jitter (carrying typical uncertainties of order 
30% - 50%) and estimated stellar rotation periods reported in 
Table |4] When both Ki » crj, and P » P^ot, conditions 
realized in the majority of the cases presented in Table |3 in- 
cluding the one studied here, little doubt can then be cast on 
the presence of the planet. If K[ > crj„ but P is similar to P^i 
(as is the case for HD 130322 and HD 192263), simultaneous 
radial-velocity, photometry, Ca II, and line bisector measure- 
ments must be carried out in order to ascertain that the planet 
is real and the observed signal is not intrinsic to the star (as 
in the case of HD 166435, also reported in Table |3 for com- 
parison). If instead P » P^ot, but Ki is comparable to cr|„ 
the situation becomes more complex, as long-term activity cy- 
cles may reproduce very low-amplitude pseudo-periodic sig- 
nals, which vary from season to season (e.g.. Gumming et al. 
1999). For example, based on the absence of a periodicity in the 
R[^^ time-series similar to that of the Keplerian signal, Hatzes 
et al. (2000) claimed HD 22049 actually harbors a planet, but 
this evidence had not been universally accepted in the past as a 
strong argument in favor of the planet hypothesis (e.g., Butler et 
al. 2003). In this specific case, the planet existence has indeed 
been proved correct (G. F. Benedict 2005, private communi- 
cation) by the recently completed HS T astrometry campaign 
(Benedict et al. 2003, 2004) on this star. 

Based on all the evidence presented here, however, the ex- 
istence of a companion around HD 8 1040 appears to be the best 
explanation for the observed radial-velocity variations. The in- 
ferred minimum mass for the object places it in the high-mass 
tail of the extrasolar planet distribution, i.e. the 18% or so of 
detected planets around nearby F-G-K stars with m2 sin ; > 5 
Mjup. This region of the planet mass distribution is populated 
by objects with minimum masses uncomfortably close to the 
widely used arbitrary dividing line between planets and brown 
dwarfs of 13 Mjup (e.g., Oppenheimer et al. 2000). As a matter 
of fact, the vanishing high-mass tail of the planet mass distribu- 
tion and the existence of a dearth of close-in (a < 5 AU) brown- 
dwarf (M < 80 Mjup) companions to solar-type stars (the so- 
called "brown dwarf desert"; e.g., Halbwachs et al. 2000) pose 
a number of puzzles to models of the formation and dynamical 
evolution of such low-mass objects, beginning with their actual 
nomenclature. Indeed, a comparatively small number of studies 
has specifically focused on explaining the existence of systems 
such as HD 168443 (the primary is orbited by two objects with 
OTsin/ = 7.7 Mjup and msini - 17.2 Mjup. Marcy et al. 2001; 
Udry et al. 2002) or HD 202206 (the primary is orbited by two 
objects with m sin / = 17.4 Mjup and m sin / = 2.4 Mjup; Correia 
et al. 2005), as opposed to the large body of work devoted to the 
study of increasingly lower-mass (M < 5 Mjup) objects, both 
single and in multi-component systems. 

High-mass planets (or low-mass brown dwarfs) such as the 
one presented here found orbiting HD 8 1040, however, are very 
interesting in their own right, as their observed properties, as 
well as the characteristics of the host stars, can provide very 
important constraints to proposed formation models. As an il- 
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lustrative example, recent models of giant planet formation by 
accretion of a rocky core (e.g., Lissauer 1993; Pollack et al. 
1996) can qualitatively reproduce the observed mass distribu- 
tion of extrasolar planets (Alibert et al. 2005; Ida & Lin 2004, 
2005), and particularly for M < 5 Mjup, down to the Neptune- 
mass (and lower) regime of some recently discovered objects 
(McArthur et al. 2004; Butler et al. 2004; Santos et al. 2004; 
Rivera et al. 2005). On the other hand, the alternative disc- 
instabiUty model (e.g., Boss 2001, 2005; Mayer et al. 2004) 
predicts that objects formed by this mechanism should prefer- 
entially populate the high-mass tail (M > 5 Mjup) of the planet 
mass distribution (Rice et al. 2003; Rafikov 2005). With im- 
proved statistics of massive planetary companions and their 
properties, as a complement to more refined theoretical stud- 
ies of the efficiency of high-mass giant planet formation (and 
of their relative frequency with respect to lower-mass planets), 
the actual roles of the two proposed formation modes could be 
better understood. 

Finally, given the projected separation of 62 mas, and given 
the relatively large mass and presumably rather young age of 
the companion, HD 8 1040 could be an interesting candidate for 
observations with future, direct near- and far-infrared imaging 
surveys of wide-separation giant planets (e.g.. Burrows 2005, 
and references therein). Also, the magnitude of the inferred as- 
trometric signature might also make this system an attractive 
target for high-precision ground-based as well as space -borne 
astrometric surveys (e.g., Sozzetti 2005, and references therein) 
which will come online in the near future. 
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